Low energy magnetic excitations of the Mni2-acetate spin cluster observed by neutron 

scattering 
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We performed high resolution diffraction and inelastic neutron scattering measurements of Mni2- 
acetate. Using a very high energy resolution, we could separate the energy levels corresponding to 
the splitting of the lowest S multiplet. Data were analyzed within a single spin model (S=10 ground 
state), using a spin Hamiltonian with parameters up to 4*'' order. The non regular spacing of the 
transition energies unambiguously shows the presence of high order terms in the anisotropy (D= 
-0.457(2) cm~^, 64= -2.33(4) 10~^cm~^). The relative intensity of the lowest energy peaks is very 
sensitive to the small transverse term, supposed to be mainly responsible for quantum tunneling. 
This allows an accurate determination of this term in zero mag netic field (B|=±3.0(5) 10"^ cm ). 
The neutron results are discussed in view of recent experiments and theories. 
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Studies of molecular nanomagnets are rapidlly grow- 
ing. The most attractive clusters consist of a few (typ- 
ically 10-20) paramagnetic ions coupled by exchange in- 
teractions, therefore at the borderline between quantum 
and classical behavior. Besides their great theoretical in- 
terest, the study of these magnetic molecules could have 
important practical consequences, since it helps to deter- 
mine the size limit for information storage. Mni2-acetate 
is the best studied spin cluster so far. In a simple ionic 
picture, an external ring of eight Mn III ions (82=2) sur- 
rounds a tetrahedron of four Mn IV ions (Si =3/2), giving 
the molecule a plate- like shape with tetragonal symmetry 
101 . Above lOK, a reversal of the molecular magnetization 
occurs by thermally overcoming the energy barrier due 
to uniaxial anisotropy (superparamagnetic behavior) . 
Below 2K, this reversal is governed by a quantum tun- 
neling mechanism as shown by a finite relaxation time in 
the magnetization at T=0 Mni2-ac also exhibits as 
a spectacular effect, regular steps in the hysteresis cycle 
Q. These observations have been interpreted by tun- 
neling processes within the ground state of the molecule 
As a peculiarity of this cluster, the tunneling pro- 
cess is thermally assisted. Mni2-ac is usually described as 
a single spin with S=10 ground state, split by anisotropy 
terms into sublevels with (-10<M<10). This picture re- 
lies on the assumption that the anisotropy energy is sig- 
nificantly smaller than the exchange energy. The S=10 
ground state corresponds to a ferrimagnetic spin arrange- 
ment of the two types of Mn ions On this basis, 
high field-high frequency EPR measurements were an- 
alyzed by considering spin Hamiltonian parameters up 
to 4*^* order ||l^. This gave the first experimental evi- 
dence of the transverse term searched for to explain the 
quantum tunneling. However, depending on the range 
of frequency or applied field considered, different sets 



of parameters were reported [^0. Very recently, the 
position of the first three excited sublevels of the S=10 
ground state was determined by optical spectroscopy in 
the far infrared (IR) range [|l2|. 

In this work, we present a detailed study of the energy 
sublevels of the ground state by inelastic neutron scatter- 
ing (INS). Neutron scattering is the most powerful tool 
to investigate spin structure and spin excitations. In con- 
trast to EPR, it is a zero field experiment, which requires 
no assumption about the Lande factor g. INS has been 
used to determine exchange and anisotropy splittings in 
numerous clusters of transition and rare earth metal ions 
p3[ . Very recently, it was used to determine the zero field 
sphtting in a Fes cluster |Q. In Mni2-ac, we previously 
studied spin excitations up to 12 meV ||l6|. We observed 
several excitations attributed either to the anisotropy or 
to the exchange terms, and determined their dynami- 
cal form factor. The position of the anisotropy level at 
low temperature (1.24 meV) agreed with the EPR de- 
termination, but we were not able to resolve the detailed 
structure of the other sublevels at higher temperature. 

We present here results obtained with a very high en- 
ergy resolution, which enable us to separate the sublevels 
of the ground state. The non regular spacing of the tran- 
sition energies unambiguously shows the presence of high 
order terms in the spin Hamiltonian. We determine these 
parameters up to 4*'' order within the single S=10 ground 
state model. Interestingly, the intensities of the lowest 
energy excitations associated with levels near the top of 
the barrier, are very sensitive to the value of the trans- 
verse 4*'' order term. This allows an accurate determi- 
nation of this very important term, although its value is 
very small. 

A [Mni20i2(CD3COO)i6(D20)4].2CD3COOD.4D20 
powder sample was synthetized using deuterated sol- 
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vents, since hydrogen has a large incoherent cross section 
which is the main background in neutron experiments. A 
deuteriation of 93% was achieved, as shown by the analy- 
sis of the neutron diffraction data. The crystal structure 
was tested between lOK and 290K by neutron diffrac- 
tion on the diffractometer G6.1 at the Laboratoire Leon 
Brillouin. Data were analyzed with the FULLPROF pro- 
gram JfTj l . The structure keeps the tetragonal symmetry 
in the whole temperature range. Above 130 K the frame- 
work of the molecule, consisting of Mn, C and O atoms 
practically coincides with that determined by single crys- 
tal X-ray diffraction . Only slight displacements of the 
acetic acid and water molecules of solvation were found. 
We precisely stated the positions of the H/D atoms, us- 
ing common constraints for interatomic distances H-0 
and H-C, and angles H-O-H and H-C-H. Moreover, we 
determined the unknown H positions in the acetic acid 
molecules of solvation. Only the water molecules contain 
hydrogen (about 33.5%), whereas the rest of the cluster 
is fully deuterated. Below 130 K, we found a slight dis- 
tortion of the cluster framework caused by displacements 
of the bridging acetate ligands. The stronger distorsion 
of the octahedral oxo coordination reduces the local sym- 
metry of Mn III with respect to room temperature, the 
environment of Mn IV being preserved (besides the usual 
lattice contraction). 

INS measurements were performed with the high en- 
ergy resolution time of flight spectrometer IN5 at ILL, 
in the temperature range 1.5K-35K. The incident wave- 
length was 5.9A. The resolution of the experiment at 
zero energy transfer (IIWIIM=27.5/iey)was determined 
using a vanadium standard. In the energy range studied 
(-1.5,1.5) meV, we checked that the q dependence of the 
neutron cross section is the same for all the inelastic tran- 
sitions. It varies very little with temperature in the range 
(0.2, 2) and agrees with previous results obtained by 
per forming constant-q scans on a triple axis spectrometer 
|l6| . We therefore added the contributions of all detec- 
tors, which strongly increased the statistical accuracy. 
The experimental spectra are shown in Fig.l for three 
temperatures. At T=1.5K, a well defined excitation is 
observed at 1.24 meV on the sample energy gain side 
{huj>0), with a peak width limited by the experimental 
resolution. It is readily attributed to an excitation from 
the lowest energy level (M=±10) to the first excited level 
(M= ±9). The peak at 0.9 meV, much less intense than 
in previous samples , was not considered in the anal- 
ysis. We also note a spurious temperature independent 
signal at the foot of the elastic peak. It was fitted by a 
gaussian and substracted. As the temperature increases, 
higher energy levels are populated, so that new excita- 
tions appear on both sides of the spectrum. At 23. 8K, 
we distinctly observe 14 well separated peaks, together 
with much smaller ones around ±0.2 meV (see the in- 
set of Fig.l). The position of the well resolved peaks is 
temperature independent. Their width, close to the res- 
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FIG. 1. Energy spectra at 3 temperatures (raw data). The 
intensity units are the same but the scales are different for all 
spectra. In inset, the low energy range on the sample energy 
loss side. 



olution limit, increases with temperature of about 25% 
between 1.5K and 25K, suggesting a finite life time. The 
energy intervals between adjacent peaks in the spectrum 
show a small but significant decrease with decreasing en- 
ergy transfer. Below 25K, the temperature evolution of 
the peak intensities reflects the depopulation of the low- 
est level at the benefit of the excited levels within the 
ground state manifold. Above 25K, all peaks broaden 
and the total inelastic intensity decreases, as the levels 
attributed to higher spin states become populated. 

In the approximation of a S=10 single spin Hamilto- 
nian, the magnetic intensity is written as fl^ ]: 
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FIG. 2. Energy sublevels of the S=10 ground state versus 
B|, as calculated for D= -0.457cm"\ B2= -2.33 10"^cm"^ 



S(Q,?ic^)-N(7^re/2)Vf^(Q)EP'.|<f|S±|i>|''5(?i^-(E/- 
E,)), with |<f|Si|i>|2=l/3(2| < f\S,\i > |2 

+ \ < f\S+\^ > \H\ < f\S^\i > \^) (1) 

is the spin component perpendicular to the scatter- 
ing vector Q. pi=exp(-Ei/kBT)/^exp(-Ei/ksT) is the 

i 

Boltzmann factor of the level i.The delta function is con- 
voluted by the experimental resolution function. The 
eigenstates |i> have the energies Ej obtained by a diag- 
onalization of the spin Hamiltonian. We adopt the spin 
Hamiltonian used to interpret EPR ]lO| , magnetization 
| p5[ and far IR experiments 

H = D[Sf_ - 1/3S{S + 1)] + B^Ol + BjOj (2) 

with O0=35St-[30S(S+l)-25]S2-6S(S+l)+3S2(S+l)2 
and Oi=l/2{SX + St). 

From equation (1), it follows that only the transitions 
with AM=0 or ±1 are allowed. When only the diag- 
onal terms are considered in equation (2) (B|=0), the 
S=10 manifold splits into 11 sublevels corresponding to 
M=±10, 9,.. O.The ±M degeneracy is not lifted, and 
when all possible levels are populated we expect nine AM 
±1 transitions on either side of the spectrum. Seven of 
these correspond to the well resolved peaks observed at 
15K and 23. 8K in Fig. 1. The decreasing energy spacing 
with decreasing energy transfer arises mainly from the 

As a first step in our quantitative data analysis in 
terms of eqs (1, 2), we neglected the transverse term B|, 
known to be small. The positions and relative intensities 
of the five most intense excitations are well reproduced, 
and the fitted values of the parameters (D= -0.457(2) 
cm~^, -2.33(4) 10~^cm~^) agree with the determi- 
nation of Barra et al by EPR |l^ and Mukhin et al by 
far IR spectroscopy p^. However, this simphfied model 
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FIG. 3. Energy spectra: filled circles are data corrected for 
the gaussian signal observed at 1.5K near zero energy transfer. 
The intensity calculation (thin line), lorentzian background 
(dashed line), elastic intensity (long-dashed line), and the sum 
of all components (thick line) are also plotted. At the top of 
the figure, intensity calculations for 3 B4 values. 



overestimates the intensities of the excitations below 0.3 
meV, which involve energy levels close to the top of the 
anisotropy barrier. 

In a second step we therefore took the off-diagonal term 
B| into account. A non zero B4 induces a mixing of en- 
ergy levels with different M that becomes more and more 
effective on approaching the top of the barrier. This is il- 
lustrated in Fig 2 where the positions of the energy levels 
Ei are plotted versus B4, for given B4 and D. The ener- 
gies of the highest excited states strongly depend on B4 
and the degeneracy of the M=±2 and ±4 levels is lifted 
for non-zero B|. At the top of Fig 3, we show the cal- 
culated intensities at 23. 8K for three different B| values, 
using the experimental energy resolution, in the energy 
transfer range (-0.5, 0.5) meV. In this range the position, 
shape, and intensities of the inelastic peaks strongly de- 
pend on B|, whereas the other peaks remain unaffected. 
Depending on the B| value (3 or 4 10~^cm~^), the energy 
transfer at 0.2 meV corresponds either to a minimum or 
to a maximum intensity, respectively. Therefore, only 
B4 values in a very narrow range can explain the exper- 
imental behavior. We find that B|=± 3.0(5) lO^^cm^^ 
gives the best agreement with the data. This value is 
slightly lower and more accurate than the one reported 
by Barra et al (B|=±4(l) lO^^cm^^). The sign of B| 
is undetermined, owing to the symmetry of the Hamilto- 
nian. In Fig 3, we present the results of the model calcu- 
lation with the best fitted values D= -0.457cm^^, B4= 
-2.33 lO-'^cm-i and B|= ±3.0(5) lO^^cm^^. In order 
to compare with experimental data at 23. 8K, we added a 
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lorentzian quasielastic background, probably due to hy- 
drogen. The agreement with experiment is excellent. We 
note a small reduction of the observed intensities for the 
excitations near 1.2 meV . It is ascribed to a partial pop- 
ulation of higher energy cluster spin states. The single 
spin S=10 model breaks down at higher temperatures 

The main result coming out of our analysis is the pre- 
cise determination in zero field of the coefficient B4 in 
the transverse term. The presence of non-diagonal terms 
in the spin Hamiltonian has been searched for by many 
theoreticians |^-|^ and experimentalists, since only com- 
ponents which do not commute with Sz could induce tun- 
nehng. In zero field, a term C(S^ -|- Si) is the lowest 
order spin Hamiltonian allowed by tetragonal symmetry. 
Such a term allows tunneling with AM=±4, as observed 
experimentally in the magnetization. Thus, our deter- 
mination of an accurate B4 value in the absence of a 
magnetic field has a direct relevance for a quantitative 
understanding of the dominant tunneling process. How- 
ever, other transverse terms are required to explain the 
tunneling transitions with AM=±1. Politi et al. have 
also argued that if only this 4*'' order transverse term 
was involved, a very small field like the earth field would 
destroy tunnelling. Other transverse terms deriving from 
dipolar coupling |2^, hyperfine interactions [|l5|, or 
a spin-phonon coupling mechanism Q , Q have been sug- 
gested in addition to the previous one. Up to now, they 
have not been observed experimentally, and our results 
provide no evidence for their existence. 

The analysis presented above is based on the assump- 
tion of a "single spin" ground state. This simple model 
describes the low energy excitations presented in this 
work very well. However, our earlier INS experiments at 
higher energies [ll6|l show that excitations involving other 
spin states are rather close to the S^IO ground state. 
Recent measurements by magnetization |l5|, or heat ca- 
pacity 1^ also point out the influence of excited states 
with different S values. Several microscopic descriptions 
of the S=10 ground state have been proposed [^ , 
1^^ , which involve different coupling schemes between 
the individual spins. Due to the large number of cluster 
spin states, 10^ in the general case, simplifications must 
be made by considering the hierarchy of the exchange 
interactions between the individual spins. Very recently, 
a new formalism for the energy has been proposed [p2[ . 
Following the " Florentine" coupling scheme M , the Mni2 
cluster is described by four dimers Mn^^^-Mn + with spin 
s=l/2, and four Mn'^+ ions (S=:2), coupled by isotropic 
exchange interactions. The dominant exchange coupling 
is between the spins within each dimer, as first sug- 
gested by j^], and in agreement with recent magnetic 
measurements in the MegaGauss range |2^. The Hamil- 
tonian includes relativistic anisotropic interactions, the 
Dzyaloshinskii-Moriya ones being the most important. 
This model qualitatively explains the spin excitations 



previously measured by INS up to 12 meV ||T^. We be- 
lieve that our new experimental results will be a crucial 
checkpoint for further theoretical work. 
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